Transverse momentum balance of dijets in Xe+ Xe collisions at the LHC 
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We present a theoretical study of the medium modifications on the pr balance (xy) of dijets 
in Xe+Xe collisions at \/snn = 5.44 TeV. The initial production of dijets is carried out by the 
POWHEG+PYTHIAS8 prescription, which matches the next-to-leading order (NLO) QCD matrix 
elements with the parton shower (PS) effect. The in-medium evolution in nucleus-nucleus collisions 
is described by the SHELL model with a transport approach. The theoretical results of the dijet 
ay in Xe+Xe collisions exhibit more imbalanced distributions than that in p+p, consistent with 
the recently reported ATLAS data. By utilizing the Interleaved Flavor Neutralisation, an infrared- 
and-collinear-safe jet flavor algorithm, to identify the flavor of the reconstructed jets, we classify 
dijets processes into three categories: gluon-gluon (gg), quark-gluon (qg) and quark-quark (qq), and 
investigate the respective medium modification patterns and fraction changes of the gg, qg, and 
qq components of the dijet sample in Xe+Xe collisions. It is shown that the gg component plays 
a key role in the increased imbalance of the dijet xj, and especially the qig2 (quark-jet-leading) 
dijets experience more significant asymmetric energy loss than the giq2 (gluon-jet-leading) dijets as 
traversing the QGP. By comparing the A(x) of inclusive, cé and bb dijets in Xe+Xe collisions, we 
observe A(r3)ina. > A(xs)ce > A(xs) pp. Moreover, pxe,pp, the ratios of nuclear modification factors 
of dijets in Xe+Xe to that in Pb+Pb, are calcualted, which indicates that the yield suppression 
of dijets in Pb+Pb is more pronounced than that in Xe+Xe due to the larger radius of the lead 


nucleus. 


I. INTRODUCTION 


The ultra-relativistic heavy-ion collisions at the Large 
Hadron Collider (LHC) and the Relativistic Heavy Ion 
Collider (RHIC) provide a unique arena to search for 
the new form of nuclear matter, the quark-gluon plasma 
(QGP), in which the degrees of freedom of the quarks 
and gluons in the protons and neutrons are released |1- 
5]. The strong interactions between the hard-scattered 
partons with the medium, referred to as the “jet quench- 
ing” phenomenon, open up new avenues to understand 
the properties of such a fantastic strongly-coupled quark 
matter [6-9] and test the fundamental theory of quan- 
tum chromodynamics (QCD) at the extremely hot and 
dense conditions [10-16]. In the past two decades, a se- 
ries of tools have been extensively investigated to reveal 
this partonic strong interaction, such as the suppression 
factor Raa of high-pr hadron/jet [17-22], the momen- 
tum asymmetry of dijets [23-34], correlations of the vec- 
tor boson associated jets (y/Z°+jets) [35-40], the global 
event geometry [41, 42] and the jet substructures [43-52]. 

Since the dijets are the dominant QCD processes in 
the hadron collisions in the experiment and have less 
influence from the underlying background, they show 
a unique glamour in jet physics. In vacuum, the par- 
ton shower effects and higher-order QCD processes may 
break the symmetry of the final-state dijets, which will 
lead to a deflection from the back-to-back azimuthal an- 
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gle and an unequal transverse momentum between the 
leading and subleading jets. In A+A collisions, because 
the two jets usually experience asymmetric energy loss 
as traversing the QGP medium, the transverse momen- 
tum balance of dijets xj = pr,2/pr,ı [53], defined as the 
ratio of the subleading to leading jet pr, can be further 
modified by the in-medium interactions and show a dif- 
ferent sensitivity to the path-length dependence of jet 
quenching [25] and jet-by-jet fluctuations of jet-medium 
interactions [54]. And more imbalanced xyz distributions 
of dijets have been observed in Pb+Pb collisions relative 
to ptp at \/snn = 2.76 TeV and \/snn = 5.02 TeV by the 
ATLAS [23, 53, 55] and CMS Collaborations [28, 56, 57], 
which are extensively investigated by the theoretical cal- 
culations [26, 27, 31-34]. 

Recently, the ATLAS Collaboration has measured the 
dijet xy in Xe+Xe collisions at \/snn = 5.44 TeV for 
the first time [58], however the timely theoretical stud- 
ies are still in lack. Because the xenon nucleus has a 
smaller radius than that of lead, studying the dijet pro- 
ductions in different collision systems will deepen our un- 
derstanding of the system size dependence of jet quench- 
ing effect [53, 55, 59-63]. Furthermore, since the dijet 
events consist of gg, qg and qq components, while the 
jet energy loss is closely related to the flavor of hard 
partons (AE,/AE, ~ Ca/Cr) [2, 18], it is of signifi- 
cance to figure out their respective modification patterns 
and assess what roles they played in the overall medium 
modifications of the dijet xy. Furthermore, the mas- 
sive heavy quarks are believed to lose less energy than 
the light quarks due to the “dead-cone” effect [64-67], 
which leads to a mass hierarchy of partonic energy loss 
AE, > AE, > AF, [49, 68-70}. It is of particular in- 


terest to explore such mass dependence of the medium 
modification on the dijet x; by the comparisons between 
the light- and heavy-flavor (such as cé@ and bb) dijets in 
high-energy nuclear collisions. 

This paper presents the first theoretical study of 
medium modifications on the dijet pr-balance xj in 
Xe+Xe collisions. The initial production of dijets is 
carried out by the POWHEG+PYTHIA8 prescription, 
which matches the next-to-leading order (NLO) QCD 
matrix elements with the parton shower (PS) effect. The 
transport approach describes dijets’ in-medium evolu- 
tion, which considers both the elastic and inelastic par- 
tonic interactions in the quark-gluon plasma (QGP). 
Firstly, we present the theoretical results of the dijet x 
in Xe+Xe collisions at ./snn = 5.44 TeV compared with 
the recently reported ATLAS measurements. Specifi- 
cally, we will discuss the flavor and mass dependence 
of the medium modification on the dijet xj. We study 
the respective medium modification patterns and frac- 
tion changes of the gg, qg, qq, as well as the qigz and 
gigz2 components in the dijet samples in both p+p and 
Xe+Xe collisions. We will also investigate the mass effect 
of the x; modifications by comparing the Axy; of inclu- 
sive, cé and bb dijets in Xe+Xe collisions. At last, we 
present the calculated results of dijets nuclear modifica- 
tion factor in Xe+Xe at \/snn = 5.44 TeV and Pb+Pb at 
/5NN = 5.02 TeV compared to the recent ATLAS data. 


Il. THEORETICAL FRAMEWORK 


In this work, we generate the next-to-leading order 
(NLO) matrix elements for QCD di-jet processes [71] in 
the framework of POWHEG-BOX-V2 [72-74] and then 
simulate the parton shower (PS) with PYTHIA 8.309 [75] 
to produce p+p events. The CT18NLO parton distribu- 
tion functions (PDF) [76] are chosen in the computation. 
Jets are reconstructed using the anti-ky clustering al- 
gorithm and radius parameter R = 0.4 as implemented 
in the FastJet package [77]. Then, the highest two pr 
jets out of the set of jets in an event are selected as the 
dijet candidate. The leading jet transverse momentum 
pr, and subleading jet transverse momentum pry are 
required to be greater than 100 GeV/c and 32 GeV/c, 
respectively. The two jets are required to be nearly back- 
to-back in azimuth with Ad = |, — ¢2| > 77/8 and to 
be in the rapidity region |y| < 2.1. If all these conditions 
are met, the desired dijet candidate is accepted. 

We calculate the normalized x; distributions in p+p 
collisions and compare that to the ATLAS data [55] as 
shown in Fig. 1. We observe that the x distributions 
calculated by POWHEG+PYTHIA8 give decent descrip- 
tions at small value of xy for all six different pp intervals 
compared to the ATLAS data, except for a little overes- 
timation of xj distributions than the ATLAS data near 
1. At each pr interval, the xy distribution peaks near 
xj ~ 1, where the leading and subleading jets are almost 
balanced. However, with the higher-order perturbative 
QCD corrections and the splitting processes during the 
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FIG. 1: (Color online) Normalized x; distributions of dijets 
in p+p collisions at ys = 5.02 TeV for six pr,ı intervals: [100, 
112], [112, 126], [126, 141], [141, 158], (158, 171] and (171, 200] 
GeV/c are compared with the ATLAS data [55]. 


parton shower in vacuum, a considerable fraction of di- 
jets are produced with imbalanced transverse momentum 
at the smaller x; region. 

The in-medium evolution of both light- and heavy- 
flavor dijets are simulated by the SHELL model [34, 78- 
85], which takes into account the elastic and inelastic 
partonic energy loss within the hot /dense QGP medium. 
Since the propagation of massive partons in the QCD 
medium can be viewed as “Brownian motion”, the trans- 
port of heavy quarks can be well described by the modi- 
fied Langevin equations, 


A(t) = PO ae (1) 
AP) = —nppAt + &(t)At — p(t). (2) 


These two equations describe heavy quarks’ position 
and momentum updates as traversing the QGP medium. 
Np is the drag coefficient controlling the energy dissi- 
pation strength of heavy quarks in the medium. The 
stochastic term &(t) denotes the random kicks as heavy 
quarks scatter with the thermal particle, which obeys a 
Gaussian distribution. The diffusion coefficient « could 
be related to the np by the fluctuation-dissipation the- 
orem K = 2npET. Note that the first two terms at the 
right-hand side of Eq. 2 represent the collisional energy 
loss of heavy quarks. In contrast, the last term —p, is 
the momentum correction caused by the medium-induced 
gluon radiation. In our framework, the Higher-Twist 
[65, 86-88] formalism has been employed to simulate the 
medium-induced gluon radiation of jet partons, in which 


the gluon radiation spectra of an energetic parton in the 
QGP can be obtained, 
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where x and k, denote the radiated gluon’s energy frac- 
tion and transverse momentum. P(x) is the QCD split- 
ting function for the splitting processes g > g +g and 
a(Q) > (Q) +9 [89], Tf = 2Be(1 — 2)/(k2 + 22M?) 
the formation time of the daughter gluon. ĝ denotes 
the general jet transport parameter in the QGP [90]. 
The last term in Eq. 3 represents the suppression fac- 
tor results from the “dead-cone” effect of heavy quarks 
(64, 91], which suppresses the probability of gluon radi- 
ation within a small cone (@ ~ Mg/E). The collisional 
energy loss is generally dominant for low-energy heavy 
quarks due to the “dead-cone” effect. In contrast, the 
radiative energy loss is usually expected to become signif- 
icant at pe > 5mg [92]. For massless partons, the colli- 
sional energy loss is estimated by the pQCD calculations 
within the Hard-Thermal Loop approximation [93, 94], 
while their radiative contribution by the same Higher- 
Twist formalism as for massive partons. The hydrody- 
namics time-space evolution of the QGP medium is de- 
scribed by the CLVisc programs [95, 96], which provides 
the temperature and velocity of the expanding hot /dense 
nuclear matter. The SHELL model has been successfully 
applied in the study of heavy-flavor jets in high-energy 
nuclear collisions, which gives satisfactory descriptions on 
a series of experiment measurements, such as pr imbal- 
ance [34], radial profiles [78-80] and fragmentation func- 
tions [84] of heavy-flavor jets, correlations of Z°+ HF 
hadron/jet [82, 83]. 


HI. NUMERICAL RESULTS AND DISCUSSIONS 


In Fig. 2(a), we firstly show the dijet xy distributions 
in Xe+Xe collisions at \/snn = 5.44 TeV calculated by 
the SHELL model as a comparison to the ATLAS data 
for four centrality bins (0 — 10%, 10% — 20%, 20% — 40% 
and 40% — 80%) and two py; intervals ([100, 126] GeV/c 
and [158, 199] GeV/c). We find that the calculations 
by the SHELL model give a satisfactory description of 
the recently reported ATLAS data for almost all four 
centrality bins and two py, intervals, only overestimate 
the data at xj ~1 for 158 < pri < 199 GeV/c in the 
central 0 — 10% collisions. Furthermore, our theoretical 
results show a more balanced x; distribution for higher 
pr dijets (158 < pri < 199 GeV/c) than the lower one, 
which is consistent with the trend observed in the AT- 
LAS measurements. To figure out the collision centrality 
and jet pr dependence of the medium modification on 
the dijet xj distributions, we also compare the dijet xj 
distributions in Xe+ Xe collisions with their p+p baseline 
for two pr intervals as shown in Fig. 2(b), and the ratios 
of XeXe/pp are shown in the bottom panels. The zxy dis- 
tributions are observed to shift towards smaller xj values 
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FIG. 2: (Color online) (a) Calculated normalized dijet £y 
distributions in Xe+Xe collisions at \/snn = 5.44 TeV are 
compared with the ATLAS data for four centrality bins (0 — 
10%, 10% — 20%, 20% — 40% and 40% — 80%) and two pr,ı 
intervals (green: [100, 126] GeV/c, red: [158, 199] GeV/c). (b) 
Comparisons of the normalized dijets x; distributions between 
Xe+Xe and p+p collisions at \/snn = 5.44 TeV for two pri 
intervals (left: [100, 126] GeV/c, right: [158, 199] GeV/c), 
and the ratios of XeXe/pp are shown in the bottom panels. 


in Xe+Xe collisions compared to their p+p baseline, es- 
pecially for the most central collisions. The phenomenon 
that the dijet transverse momentum gets more imbal- 
anced in A+A collisions due to the asymmetric energy 
loss suffered on the leading and subleading jets as the 
dijet traversing the QGP [25, 54] has been observed by 
the previous measurements in Pb+Pb collisions at 2.76 
TeV [53] and 5.02 TeV [55]. In the central collisions, such 
asymmetric energy loss between the leading and sublead- 
ing jets can be more significant due to the larger medium 
size and higher temperature than the peripheral case. In 
addition, we observe slightly weaker x; modification for 
higher pr dijets than the lower one, which is consistent 


with the expected. We can imagine that the x, of dijets 
with very high pr is relatively difficult to be influenced 
by the medium effects compared to the one with low pr. 
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FIG. 3: (Color online) Calculated normalized «x; distribu- 
tions in p+p (top) and 0 — 10% Xe+Xe (middle) collisions 
as well as their ratios (bottom) for: inclusive, gluon-gluon, 
quark-gluon, and quark-quark dijets (left column); quark-jet- 
leading and gluon-jet-leading of quark-gluon dijets (right col- 
umn). The xj distributions of quark-gluon dijets (dashed 
green line) in the right column are the same as in the left 
column. 


The flavor dependence of the medium modification of 
dijet pr balance is also an interesting topic in heavy-ion 
collisions. Due to the different color factors (C4 = 3 
for gluons, Cp = 4/3 for quarks), the gluon-initiated jets 
are expected to lose more energy than the quark-initiated 
one as traversing the hot and dense nuclear matter, and 
the “dead-cone” effect also leads to smaller energy loss 
of heavy quarks than that of the light one [2, 18, 64-67]. 
An optimal step toward this goal is determining the flavor 
of the selected jets in nucleus-nucleus collisions with an 
infrared and collinear (IRC) safe jet algorithm. In recent 
years, four IRC-safe jet-algorithm-based approaches have 
been developed to identify the flavor of the underlying 
hard partons for the final-state reconstructed jets, called 
“flavor-kp” [97], “flavor anti-kr” [98], “flavor-dressing” 
[99] and “Interleaved Flavor Neutralisation (IFN)” [100]. 
In this work, we employ the IFN algorithm to identify 
the flavor of the reconstructed jets at parton level for 
the gluon-jets, light-quark jets and the heavy-flavor jets 
in both p+p and nucleus-nucleus collisions. Note that, 
since the flavor information of the jets is accessible at 
each stage of the clustering sequence, the IFN algorithm 
provides a consistent flavor identification both to the full 
jets and their substructures. 

By utilizing the IFN algorithm, the inclusive dijet 


events can be classified into three categories: gluon-gluon 
(gg), quark-gluon (qg), and quark-quark (qq). In partic- 
ular, by distinguishing the flavor of the leading jet, the qg 
dijets can be further divided into q1g2 and giq2, denoting 
the quark-jet-leading and gluon-jet-leading quark-gluon 
dijets, respectively. In the left column of Fig. 3, we show 
the normalized x; distributions of gg, qg, qq, and inclu- 
sive dijets in both p+p and 0 — 10% Xe+Xe collisions at 
VSNN = 5.44 TeV, and the ratio of XeXe/pp is plotted 
in the bottom panel. We find that gg dijets have a more 
imbalanced initial distribution than gq. Quark and gluon 
experience different parton shower processes in vacuum 
due to their different color factors and splitting functions 
(75, 101]. In the bottom panel, it is observed that the qq 
and gg dijets have almost the same medium modifica- 
tion, while the qg dijets have the stronger suppression 
near xj ~ 1 and enhancement near x; ~ 0 than the oth- 
ers. The different color charge carried by the leading and 
subleading jets of qg dijets leads to enhanced asymmet- 
ric energy loss in Xe+Xe collisions relative to the qq and 
gg one. However, we note that the qg dijets contain two 
types of subsets, q1g2 and gıq2, which may have differ- 
ent modification patterns. In the right column of Fig. 3, 
we observe stronger enhancement at small xj and also 
stronger suppression at xj ~1 on the x; distribution of 
Mg2 dijets than that of gıq2. Compared to the gq di- 
jets, the leading jet of q1g2 dijets loses less energy while 
the subleading one loses more. In other words, the flavor 
configuration of q1g2 dijets makes them experience more 
significant asymmetric energy loss than g;q2 as traversing 
the QGP medium. 


On the other hand, we also estimate the component 
fractions of dijet sample in p+p and 0— 10% Xe+Xe col- 
lisions at \/snn = 5.44 TeV shown in Fig. 4, which would 
be helpful to understand the role of jet flavor played 
in the jet-medium interactions. In the left column, we 
show the fractions of gg, qg, and qq in the inclusive di- 
jets in p+p (top) and Xe+Xe (middle) collisions, as well 
as their differences (bottom). Firstly, we find that the qg 
has the most prominent initial fraction (~ 50%) in p+p 
collisions, and the fraction is increased at small xj but 
reduced at xj ~ 1 in Xe+Xe collisions. Second, the frac- 
tion of gg is overall reduced. At the same time, that of 
qq is enhanced because, generally, the gg dijets lose more 
energy, which makes it more difficult to survive in jet 
selection relative to qq. Since the gg has the most signif- 
icant fraction in the dijet sample, the enhanced fractions 
of qg at small x; in Xe+Xe is the key point that leads 
to the increased pr imbalance of inclusive dijets. Fur- 
thermore, it is essential to address the fraction changes 
of the qıg2 and gıq2 subsets in A+A collisions. In the 
right column of Fig. 4, it is observed that the fractions 
of q1g2 and gıq2 have opposite behaviors in Xe+Xe col- 
lisions compared to their initial values. The fraction of 
M92 is significantly enhanced at small x; after traversing 
the QGP, while that of gig is decreased at this region. 

To quantify the overall shift of the xj distribution in 
Xe+Xe collisions relative to p+p, we calculate the aver- 
age values ((xy)) of dijet x3 distributions and their differ- 
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FIG. 4: (Color online) Calculated fraction of subset dijets in 
inclusive dijets as the function of x; in p+p (top) and 0— 10% 
Xe+Xe (middle) collisions as well as their differences (bot- 
tom) for: gluon-gluon, quark-gluon, and quark-quark dijets 
(left column); quark-jet-leading and gluon-jet-leading quark- 
gluon dijets (right column). The fraction distributions of 
quark-gluon dijets in the right column (dashed green line) 
are the same as in the left column. 


ences (A(x3)) between p+p and A+A collisions, defined 


as follows 
= 1 dNpair 
(x3) = J Noan day xydxy, (4) 
A(x) = (&3)pp — (TI) Aa. (5) 


To address the mass dependence of the medium mod- 
ification of dijet xj, in the left column of Fig. 5 we 
show the (x) of inclusive, cé and bb dijets in p+p and 
Xe+Xe collisions as a function of centrality, as well as 
their differences A(xy). We find that A(z) decreases 
monotonously from central to peripheral collisions, as 
we found in Fig. 2(b). It is also observed that the 
A(ay) of these three kinds of dijets obey the hierarchy 
A(xj)ina. > A(x3)ce > A(xy)p5 in Xe+Xe collisions for 
the same centrality bin. It indicates that the massive di- 
jets suffer less asymmetric energy loss in A+A collisions 
than the massless light flavor one. The future measure- 
ments focusing on these comparisons will help test the 
mass effect of jet energy loss. In the right column of 
Fig. 5, A(ay) of qıg2 and giq2 dijets in Xe+ Xe are also 
plotted, and the former has significantly large values than 
the latter for each centrality bin. 

To quantitatively characterize the relative dijet yield 
suppression between Xe+Xe and Pb+Pb collisions, the 
nuclear-modification factors for leading jets, pxe,pp(pr.1) 
is defined as follows (similarly, pxe,pp(pr,2) can be de- 
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dijets (left) and qig2, giq2, qg dijets (right) in p+p (top) 
and Xe+Xe (middle) collisions as a function of the collision 
centrality, as well as their differences A(x3) (bottom). 


fined for subleading jets) 


1 do’ /dpr, i 
(Taa) /dpr, 2 
Rxexe(Pr,1) (7) 
Rpppp(pr,1) | 


Raa(pr,1) = 


PxXe,Pb(PT,1) = 


In Fig. 6, we show the calculated pxe,pp(pr.1) (left) 
and pxe,pp(pr,2) (right) of dijets for three centrality bins: 
0 — 10%, 10% — 20% and 20% — 40%. Note that we 
use a cut xj > 0.32 in the calculation to be consistent 
with ATLAS’s treatment in the measurements. We ob- 
serve that the values of pxe,pp are generally larger than 
one for both leading and subleading jets for all centrality 
bins. It indicates a weaker yield suppression of dijets in 
Xe+Xe collisions than that in Pb+Pb for the same cen- 
trality bin, and this phenomenon is still evident even at 
peripheral collisions. These findings are consistent with 
the previous phenomenological studies [61-63]. Since the 
nucleus of xenon has a smaller radius than that of lead, 
the system size and the mean temperature of the QGP 
medium formed in Xe+Xe collisions is expected to be 
smaller than that in Pb+Pb collisions within the same 
centrality interval [85]. Hence, dijets traverse a longer 
path-length medium and experience more effective en- 
ergy loss in Pb+Pb collisions. 
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FIG. 6: (Color online) Calculated ratios of Xe+Xe and 
Pb+Pb pair nuclear modification factors, pxe,pp, evaluated 
as a function of pr, (left) and pr,2 (right) in the same cen- 
trality intervals, and compared to the ATLAS data. 


IV. CONCLUSION 


In this paper, we present the first investigation on 
the medium modifications of dijet pr balance (x3) 
in Xe+Xe collisions at \/snn = 5.44 TeV. The ini- 
tial xy distributions of dijets are calculated by the 
POWHEG+PYTHIA8 prescription, which matches the 
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